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We derive the exact equality, referred to as the fluctuation relation for heat engines (FRHE), that relates 
statistics of heat extracted from one of the two heat baths and the work per one cycle of a heat engine operation. 
Carnot's inequality of classical thermodynamics follows as a direct consequence of the FRHE. 

PACS numbers: 05.60.-k, 05.40.-a, 82.37.-j, 82.20.-w 



The central thermodynamic result for the efficiency of heat 
engines, which was established by Carnot in the 19th century, 
reads 
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where W is the work performed by a heat engine, Qh is the 
heat that was absorbed from the hotter bath during the cyclic 
operation, and T c and Th are the absolute temperatures, re- 
spectively, of the cooler and hotter baths, which are assumed 
constant during the process. The motive power of a heat en- 
gine originates from the passage of heat from a hotter to a 
cooler heat bath, as shown in Fig. [JJ Part of this heat flux can 
be converted to useful work by means of a cyclic operation of 
a central system that is coupled to both of the baths. When 
work is extracted in this way, there is always a heat that is 
lost in the cooler bath. Carnot's inequality sets a limit on the 
useful work that can be extracted from the absorbed power. 

The inequality (JTJ was originally derived to describe oper- 
ations of macroscopic machines. In modern times, a lot of at- 
tention has been devoted to nonequilibrium thermodynamics 
of systems so small that thermal fluctuations beyond the realm 
of the Gaussian approximation can be observed and charac- 
terized. Stochastic behavior of mesoscopic heat engines is 
a growing field of modern research [1| with applications to 
nanostructures, including nanocoolers |2| and molecular mo- 
tors (3J. Surprisingly, the stochastic behavior of strongly 
driven small systems was found to satisfy a universal con- 
straint, as articulated in the Fluctuation Theorems (FTs) |0}- 

FTs transform classical thermodynamic inequalities into 
equalities for exponents of thermodynamic variables. Eq. ([TJ 
is one of the classical inequalities. Consequently, it can be 
derived starting from a variety of known fluctuation relations 
ll6l [T0l but, to the best of author's knowledge, there is no writ- 
ten FT that relates exactly the same variables as in ([TJ in a 
single exact expression. The lack of such an equality may 
leave an impression that different formulations of the 2nd law 
of thermodynamics, one of which is based on ([TJ, are not com- 
pletely equivalent when they are considered from the point of 
view of fluctuation relations. In this article, we resolve this 
problem by deriving a FT, which we will call the FRHE, that 
relates the same variables as appear in ([TJ but now in the form 
of an exact equality. 

The FRHE states that for an arbitrary heat engine the work 
per cycle of the driving protocol and the heat absorbed from 




FIG. 1. Heat engine. Hot and cool baths (boxes) at temperatures, 
respectively, Th and T c , are coupled to the central system (circle). 
Qh is the heat that is absorbed from the hot bath and Q c is the heat 
that is released to the cool bath during the engine operation. Periodic 
changes of parameters of the central system coupled to baths and 
external fields lead to production of work, W, at a cost of heat fluxes 
from the hot to the cool baths. 



the hotter bath per cycle are related by 
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where the averaging is over many realizations of the thermo- 
dynamic cycle. Applying Jensen's inequality, exp((x)) < 
(exp(z)), Eq. (TfJ leads to (W)/T c < (Q h )(l/T c - 1/T h ), 
which is equivalent to ([TJ. Thus the primary importance of 
the FRHE is that it relates the same variables as in the Carnot's 
theorem, thus promoting the latter to the exact equality. 

Here we should specify the initial and the final conditions 
more clearly. We will assume that the central system is in 
the thermodynamic equilibrium with the cooler bath at the be- 
ginning of the protocol. The central system switches on the 
coupling to the hotter bath and performs arbitrary changes of 
parameters until finally it returns all parameters to their initial 
values, including switching off the coupling to the hotter bath, 
which ends the thermodynamic cycle. Eq. |2j becomes exact 
only for such protocols. 

There are two limiting cases for which the FRHE (|2j re- 
duces to the already known FTs. First, imagine that the system 
is completely decoupled from the hot bath. Then Q n = 0, and 
(|2j becomes equivalent to the Bochkov-Kuzovlyov-Jarzynski 
equality HUD, (exp(W/T)) = 1, for a periodically driven 
system coupled to a single heat bath. On the other hand, con- 
sider a system that does not perform any work. This happens 
e.g. when we simply couple the central system to the heat 



2 



baths without connecting the central system to other external 
fields. Then the central system works just as a heat conduc- 
tor between two baths and Eq. Q becomes equivalent to the 
so-called Exchange Fluctuation Theorem ifTTl . 

We will not discuss the general domain of applicability 
of the FRHE and restrict its derivation only to the domain 
of stochastic thermodynamics lfT2l - [2"2"l . in which continuous 
Markov chains are used to model fluctuations of physical sys- 
tems. We note, however, that there are closely related expres- 
sions, such as Eq. (9) in Ref. ifTOll that were derived starting 
from Hamiltonian equations of motion. This indicates that 
extensions of the domain of the FRHE applicability should be 
possible. 

In our derivation of the FRHE, we assume that the dynam- 
ics of the central system is described by stochastic motion 
among a finite number N of discrete states at chemical poten- 
tials Ei, i = 1, . . . , N. We assume that these discrete states 
are separated by barriers and that each barrier is strongly cou- 
pled to one of the heat baths. Each time the central system 
makes a transition from a state i into a state j via a barrier 
it exchanges energy with a corresponding bath, which forces 
the kinetic rates of transitions to have the Arrhenius form. 
This guarantees that, at thermodynamic equilibrium, the sys- 
tem state probabilities are described by the Boltzmann-Gibbs 
distribution. For a transition through a barrier, which is cou- 

h/c 

pled to the hot/cool bath, such a kinetic rate, k^ , is given 
by fck/ c — g i ^ e Ei/T h/c ^ whgj-g — gj. qj-q parameters that 
characterize barrier sizes [23 1. Every time such a transition 
occures, heat is transferred between the central system and 
a corresponding bath. The corresponding change of the en- 
ergy of the central system, Ej — Ei, is attributed to absorp- 
tion/release of energy from/to a hotter/cooler heat bath, i.e., 
respectively, either 5Qh — Ej — Ei or SQ C = Ei — Ej. 

We consider protocols such that parameters <7y and Ei are 
the same at the beginning and at the end of the parameter driv- 
ing. The probability, P F (X), of a trajectory, X, and the prob- 
ability, Pb(X), of its time-reversed counterpart, X, in the 
same system driven by a time-reversed protocol, satisfy the 
relation Q [H ED 
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Here S(X) = — lnf} r (kij(t r )/kji(t r )), where the product 
is over all transition time moments, enumerated by index r, 
along a trajectory X, and where t r is the time of the transi- 
tion from the state i into the state j along X; kij(t r ) is the 
kinetic rate of this transition and kji(t r ) is the kinetic rate of 
the transition in the opposite direction at time t r . f^ q and /° q 
are the equilibrium probabilities of having, respectively, ini- 
tial and final (at time r of the end of the protocol) states of the 
trajectory when the central system is decoupled from the hot 
bath. Reformulating ([3jl in terms of Ei and <?y, we find that 
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ferred from/to hot/cool baths to/from the central system, as- 
suming that at the beginning of the cycle Qy c = 0. 

The energy balance between the work, W, the heat fluxes, 
Qh, Qc, and the change of the internal energy of the central 
system, AE = E(t) — E(0), is generally given by 



Qh - Qc = W + AE. 



(5) 



Our goal is to relate Qh and W. One can think that since 
Eq.(|5jl is the only exact constraint among 4 variables Qh, Q c , 
W and AE, further fluctuation relations should involve at 
least 3 of those variables. However, when averages of expo- 
nents are considered, additional cancellations happen. Elimi- 
nating Q c from |5jl and substituting the result in Q we find 
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where Q 
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Sr^Qh/c(*r) are total heat amounts trans- 



where we used the fact that due to the detailed balance condi- 
tions, e £ (°)/ T ^ Cq = /eq e ^(r)/T c> 

Let p(Qh, W) be the probability density to observe given 
values of Qh and W. The standard manipulations lfT2l then 
give us 

P (Q h , W) = J dXP F (X)S(W - W F (X))S(Q h - Qi{X)) 
e-QbWT^-i/Tj-w/T. j d XP B (X)6{W + W B (X))6(Q h + 

where Q^(X) = —Q^{X) are the quantity of heat absorbed 
from the hotter bath during forward and backward protocols in 
motion along, respectively, the forward and the backward tra- 
jectories. Analogously we denote by W F {X) = — W B {X) 
the amount of work done. Summing in (|7]i over all possi- 
ble Qh and W, we obtain the FRHE which completes its 
derivation. 

In conclusion, we derived a fluctuation relation that couples 
the work performed by a heat engine with the heat flux from 
the hotter heat bath. The equality describes performance of an 
arbitrary heat engine that can be driven by arbitrarily fast and 
strong time-dependent fields. We also allowed both heat baths 
to be coupled strongly to the central system at the same time, 
except the beginning and the end of the driving protocol. The 
FRHE directly leads to Carnot's inequality for macroscopic 
heat engines and, in this sense, it may be regarded as gener- 
alization of this crucial thermodynamic principle. Certainly, 
at present, Carnot's inequality remains the true fundamental 
law of physics because the full domain of applicability of the 
FRHE in the real world remains to be understood. FTs have 
already been used for practical purposes, including the devel- 
opment of accelerated free energy sampling algorithms 11241 . 
The FRHE should similarly help the research on optimiza- 
tion and control of nanoscale device structures such as electric 
coolers and molecular motors Q. 
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